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I briefly review lattice QCD calculations that study the 0^"*" glueball and discuss implications for light flavour 
singlet 0"*"^ mesons. 



1. Introduction 

The 0"^^ glueball in pure gauge theory is a well 
defined quantity that can be used to test our un- 
derstanding and ability to solve theory. It is also 
interesting to understand how, or if, the 0^"*" glue- 
ball contributes to the physical light flavour sin- 
glet O"'"'^ mesons. 

Finding glueball degrees of freedom in the 
flavour singlet mesons is complicated, be- 
cause other non-perturbative objects such as 
tetraquark, meson molecule or, even quark- 
antiquark degrees of freedom can also be building 
blocks of scalar mesons. There are reviews that 
discuss these broader issues in more detail [1-3]. 

1.1. Background to lattice QCD 

Lattice QCD is based on a Monte Carlo process 
where a statistical sample of vacuum gauge fields 
is produced. On each sample of the QCD vac- 
uum, an interpolating operator creates a hadron 
and after a specific time interval the hadron is 
destroyed. The choice of interpolating operator 
is particularly important for hadrons where it is 
not clear how the hadron is built out of quarks 
and gluons. 

For example, to create a light flavour singlet 
O"'"'*' hadron, possible interpolating operators are 

01 = qq (1) 

02 = qibqqibq (2) 

03 = Uplaq (3) 

where Upiaq is a spatial plaquette of gauge fields 



with 0^+ symmetry, and g is a light quark oper- 
ator. 

The majority of recent lattice QCD calcula- 
tions that include the dynamics of sea quarks 
have pion masses as low as 300 MeV [2], with 
a range of volumes and lattice spacings. These 
parameters have allowed lattice QCD to make 
contact with chiral perturbation theory for light 
pseudoscalar mesons [4]. For some quantities, 
there are now precision results from lattice QCD. 
For example the HPQCD collaboration obtained 
iricimc) = 1.268(9) GeV for the mass of the 
charm quark [5]. The LHPC collaboration [6] 
computed the nucleon axial charge to be gA = 
1.212(84) from unquenched lattice QCD. The low 
lying glueball spectrum in pure SU (3) gauge the- 
ory was accurately computed nearly ten years 
ago [7,8] 

Unfortunately, lattice QCD calculations of 
scalar mesons are not as accurate as those of other 
quantities. The lattice QCD correlators for scalar 
mesons are more noisy than for p and tt mesons, 
so much higher statistics are required. The light 
scalar mesons decay via S-wave decays, and cur- 
rent lattice QCD calculations are in the quark 
mass regime where some decay channels to two 
mesons are open. The results I will present for 
flavour singlet 0++ mesons largely use the last 
generation of lattice QCD calculations that are 
quenched, or dynamical QCD calculations with 
pion masses above 500 MeV [9, 10]. 
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Method 


"^2++/^7^o++ 


lattice SU{3) [7] 
lattice SUiN) [12] 


1.39(4) 
1.46(11) 


strong coupling lattice [13] 
black hole [14] 
Klebanov-Strassler [15] 
Frasca, strong coupling [16] 


1.2 - 1.25 
1.73 
1.37 
1.33 



Table 1 

Glueball ratios from different calculations 



1.2. Glueballs as a theoretical testing 
ground 

The masses of glueballs from quenched QCD 
are a good test of calculational techniques, such 
as methods based on a gauge/gravity duality 
(see [11] for a review). Many of these new meth- 
ods produce results in the large N (number of 
colors) limit. Teper and collaborators [12] have 
studied glueballs for different SU{N) groups on 
the lattice. In general they find the correc- 
tions to the glueball masses in the large N limit 
to be small. For example [12] the mass of the 
mo++ glueball in units of string tension a was 
found to be 



^3.37(15) ■ '-''^''^ 



(4) 



Some of the methods based on gauge/gravity 
duality produce results in the strong coupling 
limit, so it is (perhaps) interesting to compare 
them with older lattice results computed using a 
strong coupling expansion in p-. It is not always 
easy to get a physical scale from all calculations, 
so I compare the dimensionless ratio of the mass 
of the 0++ and 2++ glueballs in table[Il The black 
hole and Klebanov-Strassler results are based on 
calculations using gauge/gravity that use differ- 
ent metrics. Considering that the method that 
uses the Klebanov-Strassler [15] metric has a 
massless 0+"'" state, and that the calculation by 
Frasca [16] in pure Yang-Mills has an additional 
0"^"*" state around 530 MeV, the result from the 
strong coupling lattice QCD calculation is in fair 
agreement with the SU(3) result. Brower et 
al. [14] compare in detail their results from a 
gauge/gravity duality calculation with those from 
lattice QCD [7]. 



1.3. Flavour singlet 0++ mesons 

In nature and unquenched lattice QCD cal- 
culations glueball and qq operators will mix, 
so glueballs do not exist as separate particles. 
The lightest flavor singlet 0"'""'" mesons listed in 
the PDG [17] are: /o(600), /o(980), /o(1370), 
/o(1500), and /o(1710), so it is expected that 
0"*"+ glueball degrees of freedom will contribute 
to some of these mesons. There are claims that 
the /o(980) and /o(600) may be molecules or 
tetraquark [17], so may not couple to qq inter- 
polating operators in lattice QCD calculations. 

Morningstar and Peardon [7] obtained Mq++ 
^ 1730(50) (80) MeV for the mass of the hght- 
est O"^"*" glueball from quenched QCD. Chen et 
al. [8] recently found Mo++ = 1710(50) (80) MeV 
from quenched QCD. The quark model predicts 
that there should only be two O"*""*" mesons be- 
tween 1300 and 1800 MeV, so if the mixing be- 
tween the glueball and qq operators is weak, then 
the 0++ glueball is hidden inside the /o(1370), 
/o(1500) and /o(1710) mesons. Klempt and Za- 
itsev have recently argued that the experimental 
data don't support that the /o(1370) is a separate 
resonance [18]. 

Weingarten and Lee [19] used quenched lattice 
QCD to estimate the mixing matrix between the 
glue and qq states. Weingarten and Lee [19] pre- 
dicted that the /o(1710) meson was 74(10)% 0++ 
glueball, and hence the mixing between the O"*""*" 
glueball and qq states is weak. 

There are claims [20] that continuum phe- 
nomenology is more consistent with a sizable con- 
tributions from the O"*""*" glueball to the /o(600) 
and /o(980) mesons, so it is important to study 
the effect of sea quarks on the masses of the state 
created by glueball interpolating operators. 

The SESAM collaboration studied the glueball 
spectrum on unquenched lattices [21]. McNeile 
and Michael studied the light 0+"'" spectrum with 
unquenched QCD [22] at a coarse lattice spacing 
and found the mass of the lightest flavour singlet 
0++ meson was very light. Using O"*"*" glueball 
operators. Hart and Teper [23] found that 

UN quenched = 0.85(3) Afo+ + Q„e„c;ied (5) 

at a flxcd lattice spacing of 0.1 fm. The UKQCD 
collaboration [24] separately studied 0++ glue- 
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ball and 0+"'" qq operators on improved staggered 
gauge configurations, however higher statistics 
and an analysis similar to the one by Bernard 
et al. [25] is required. 

Unfortunately, the existing unquenched lat- 
tice QCD calculations of the flavour singlet 0++ 
mesons don't have the range of lattice spac- 
ings where a continuum extrapolation can be at- 
tempted. In quenched QCD it was found that the 
lattice spacing dependence of the mass of the 0++ 
glueball was strong. The use of a Symanzik im- 
proved gauge action by Chen et al. [8] and, Morn- 
ingstar and Peardon [7], produced a smaller de- 
pendence on the lattice spacing of the scalar 0'*'"'' 
glueball mass, than for calculations that used the 
Wilson plaquette action. This is relevant to un- 
quenched calculations, because any suppression 
of the mass of the flavour singlet 0"''^ mass may 
be due to lattice spacing effects. 

The SCALAR coUaboration [26], used un- 
quenched lattice QCD, with Wilson fermions 
and the Wilson gauge action, to study the 0++ 
mesons. At a single lattice spacing a ~ 0.2 fm, 
with qq interpolating operators only, they obtain 
'n^qq fTT'p- The lattice spacing dependence of 
this result needs to quantified. 

In unquenched QCD, both glue and qq states 
will couple to singlet 0++ mesons, so it is bet- 
ter to do a variational fit with both types of 
operators as basis interpolating operators. The 
variational technique analysis of the singlet 0++ 
mesons was done by Hart et al. [27]. A combined 
fit to 0++ glue and qq interpolating operators 
with two types of spatial smearing sources was 
done. The calculation used the non-perturbative 
improved clover action at a single lattice spac- 
ing [9]. Configurations from CP-PACS [10] with 
the Iwasaki gauge action and tadpole improved 
clover action were also used in the analysis, be- 
cause this calculation should be less affected by 
lattice artifacts. A summary plot of the results, 
in units of rg (I/t'o ^ 400 MeV) is in figure [T] 
(updated from [27]). The data with the bursts 
and squares (with the pion masses written near 
them) in figure[T]shows an additional reduction of 
the mass of the 0++ state over the pure glueball 
operators, as used by Hart and Teper [23]. 

Mathur at al. [28] recently claimed to get a re- 
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Figure 1. Summary of unquenched results for 
lightest fiavour singlet 0"'""'' mesons from [27]. 
The unquenched results are from SESAM [21], 
UKQCD-I [22], and UKQCD-II [27]. 



suit for the mass of the /o(600) (cr) from quenched 
lattice QCD with pion masses as low as 180 
MeV. Using the interpolating operator in equa- 



tion [21 they obtain rrifg 



/o(600) 



550 MeV. Mathur 



et al.'s [29] calculation is discussed in shghtly 
more detail in [2]. The effect of sea quarks on 
this calculation needs to be quantified. 

In [27] an attempt was made to compute the 
decay width for /o decay to two pions. Unfortu- 
nately much higher statistics will be required to 
obtain an accurate value for that width. Recent 
unquenched lattice QCD calculations have light 
enough quarks that the two meson decays of some 
scalar mesons are allowed and some preliminary 
evidence for 0++ decay has been presented [30]. 

Perhaps a more mundane issue with the im- 
proving unquenched lattice QCD calculations of 
flavour singlet quantities is just increasing the 
statistics in the Monte Carlo estimate. In ta- 
ble [21 I show the number of measurements done 
in some quenched and unquenched calculations 
of 0+"*" glueballs. This type of comparison be- 
tween number of estimates, can be misleading 
because it depends on the autocorrelation times. 
The qualitative message from table [21 is that un- 
quenched calculation of the flavour singlet 0++ 
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Group Method statistics 

Morningstar ct al. [7] quenched 6360 

Chen et al. [8] quenched 10000 

Hart et al. [27] unquenched 500 



8. Y. Chen et al., Phys. Rev. D73 (2006) 014516, 
hep-lat/0510074, 

9. UKQCD, C.R. Allton et al., Ph ys. Rev. D65 
(2002) 054502, hep-lat/0107021, 

10. CP-PACS, A. Ah Khan et al., Ph ys. Rev. 
D65 (2002) 054505, |hep-lat/0105015l 

11. M. Caselle, Int. J. Mod. Phys. A15 (2000) 
3901, hep-th/0003 119, 

12. B. Lucini'and M. Teper, JHEP 06 (2001) 050, 
"hep-lat/0103027; 

13. J. Smit, Nucl. Phys. B206 (1982) 309, 

14. R.C. Brower, S.D. Mathur and C.I. 
Tan, Nucl. Phys. B587 (2000) 249, 
lhep-th /0003115 

15. X. A mador and E. C aceres, JHEP 11 (2004) 
022, hep-th/0402061', 

16. M. Frasca, (2008), 0807.4299, 

17. Particle Data Group, W.M. Yao et al., J. 
Phys. G33 (2006) 1, 

18. E. Klempt and A. Zaitsev, Phys. Rept. 454 
(2007) 1, 0708.4016, 

19. W.J. Lee and D. Weingarten, Phys. Rev. D61 
(2000) 014015, hep-lat/99 10008, 

20. G. Mennessier et al., C2007V larXiv:0707.4511l 
[hep-ph], 

21. TXL, G.S. Bah et al., Phys. Rev. D62 (2000) 
054503, hep-lat/0003012, 

22. UKQCD, C. McNeile and C. Michael, Phys. 
Rev. D63 (2001) 114503, hep-lat/0010019, 

23. UKQCD, A. Hart and M. Teper, Phys. Rev. 
D65 (2002) 034502, hep-lat/0108022, 

24. E.B. Gre gory et al ., PoS LAT2005 (2006) 
027, |hep-rat/05r006~6i 

25. C. Bernard et al., (2007), larXiv:0707.2402l 
[hep-lat] , 

26. SCALAR, T. Kunihiro et al., Phys. Rev. D70 
(2004) 034504, hep-ph/0310312, 

27. UKQCD, A. Hart et al., Phys. Rev. D74 
(2006) 114504, |hep-rat/060 8026, 

28. N. Mathur et al., Phys. Rev. D76 (2007) 
114505, hep-ph/0 607110 , 

29. N. Mathu r et al., PI^ s. Rev. D70 (2004) 
074508, hep-ph/0406196', 

30. ETM, .C. Michael and C. Urbach, (2007), 
lar2£iy:0709.4564 [hep-lat], 

31. E.B. Gregory et al., PoS LAT2007 (2007) 
099, 0710.1725, 

32. H.B. Meyer, (2008), 0808.3151, 



Table 2 

Comparison of statistics between quenched and 
unquenched glueball calculations. 



mesons need at least 10 times as much statistics 
as currently used. Lattice QCD calculations with 
higher statistics are definitely required to study 
unquenching on the mass obtained from 2++ glue- 
ball interpolating operator [23]. There is a high 
statistics unquenched lattice QCD calculation in 
progress of light flavour singlet mesons, that has 
generated 6000 configurations at a single value of 
the lattice spacing [24,31]. 

Chen et al. [8] and Meyer [32] have recently 
computed the matrix element of the energy mo- 
mentum tensor with glueball states in quenched 
QCD calculations. Meyer [32] uses the computed 
matrix elements to study J/?/; radiative decay. 

2. Conclusions 

There is "some" evidence that the fiavour sin- 
glet 0"'""'' qq and glueball interpolating operators, 
in unquenched lattice QCD calculations, are cou- 
pling to states around or below 1 GeV [27]. Al- 
though a continuum extrapolation is required for 
definite results. Unquenched lattice QCD calcu- 
lations with 0"*"+ tetraquark interpolating opera- 
tors are required to clarify the composition of the 
lightest 0"'"'" resonance. 
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